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ABSTRACT 


We conduct a spectroscopic search of quasars observed by the Sloan Digital 
Sky Survey (SDSS) with broad absorption line (BAL) troughs due to Mgll and 
troughs due to Fe ll that simultaneously exhibit strong Balmer narrow emission 
lines. We hud that in a redshift range of 0.4 < z < 0.9 approximately 23 of the 70 
Mg II BALs and 4 of a subset of 15 Fe ll BALs exhibit strong Balmer emission. We 
also hnd signihcant fractions of Mgll BALs (approx. 23%) and those Mgll BALs 
with Fell troughs (approx. 27%) have strong continuum reddening, E(B—V) > 
0.1. From measurements of the Balmer decrement in three objects, we hnd 
similarly signihcant reddening of the narrow emission line region in 3 of the 4 
objects; the narrow emission lines in the fourth object are not measurable. We 
also include one object in this study not taken from the SDSS sample that shows 
Fe II absorption and strong narrow emission, but due to measurement uncertainty 
and low continuum reddening the comparison is consistent but inconclusive. We 
hnd a trend in both the Mg ll and Fe ll BAL samples between the narrow emission 
line reddening and continuum reddening. Because the narrow line reddening is 
consistent with the continuum reddening in every object in the two SDSS samples, 
it suggests that the reddening sources in these objects likely exist at larger radial 
distances than the narrow line regions from the central nucleus. 


Subject headings: quasars: absorption lines, galaxies: evolution 
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Introduction 


Mass outflows in quasars, seen as blueshifted broad absorpti on lines (BALsb have 


recen tly become a plausible explanation for quasar feedback (e.g., 


Faucher-Giguere et al. 


2 OI 2 I). Feedback has been invoked to explain several phenomena such as the M—a 


relationship for supermassiv e black holes (SMBH) and their host galaxies (e.g., 


1998; 

Di Matteo et al. 

200 

5) as well as 

die Qi 

(e.g., 

Hopkins et al. 

2013; 

Dubois et al. 

2013) 


Silk fc Ree, 


uenching of star formation in the host galaxy 


201311 . The recent advancement is due in large part 


to studies of a subclass of BAL quasars that exhibit absorption troughs from both resonance 


and metastable state lines of 
FeLoBALs; 


Korista et al. 


2008 


ow ionization species, namely from Fell and Fell* (i.e ., 


Moe et al. 


2009 


Dunn et ah 


in one case from high ionization lines (S IV and S VI*; 


2010 


Bautista et al 


Borguet et al. 


2 OI 0 I) and 


2 OI 3 I) . The combination 


of these lines with photoionization modeling yields the density and therefore the distance 


and energy output of the outflows. These observational results have 


explain the p 

lysica 

Debuhr et al. 

2010; 


mechanisms 


Dehind the driving of the gas (e.g.. 


Silk fc Nusser 


ed to recent work to 


Ostriker et al. 


2010 


2010h 


Unlike the more common high ionization BALs, which dem onstrate strong abso rption 


due to C IV and app ear in approximately 20% of all quasars flKnigge et al. 


2008 


Hewett fc Foltz 


20031) ■ low ionization absorption troughs due to species like Mg ll (LoBALs) 


2000; 

Hall et al. 2002: 

Silk & Nusser 

2 OIOII. Furthermore, many 

joBALs ( 

Zhang et al 

r 

2010) 

and FeLoBALs ( 

Weymann et al. 

1991: Snravberrv & Foltz 

1992; Reichard et al. 

2003; 

Urrutia et al. 200 

8) display signs of large color excess due to dust extinction. The 


extinction commonly observed in these obj ects appears to be best described by a S mall 


Magellanic Cloud (SMC) extinction curve (IRichards et al. 


(Richards et al. 

2003; 

Hopkins et al. 

2004) 


combination of rarity and reddening properties of these objects has led studies to postulate 
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about their relevance in quasar evolution. One plausible scenario is that FeLoBALs are 
normal BAL quasars with a specialized sight-line that grazes the edge of the putative dusty 
torus (see Hall et ah 2002). This simultaneously explains the strong reddening present in 
many objects and provides a source of opacity that decreases the ionizing flux impinging 
on the outflowing gas allowing for clouds to form singly ionized species (see Dunn et ah 
2010 and references therein). An alternative sce nario suggests that FeLoBALs represe nt 


an evolutionary sta. t e during a quasars 


Hopkins et al. 


2008 


Urrutia et al. 


2008 


ifetime flEgami et a. 


Glikman et al. 


1996 


Becker et al 


19971: 


2 OI 2 I) where the quasar sheds its 


obscuring dust and gas. The dust source during this process could provide reddening at 
any distance from the nucleus. 


While recent studies of FeLoBALs have determined distances for these outflows, the 
radial location of the dust source is unknown. This was especially relevant in the object 
QSO0318—0600, where we accurately determined the density of the outflowing gas due 
to the presence of troughs from metastable, excited-state lines (Fell, Sill, and Nill). 
Combining the density with photoionization modeling yielded a radial distance of several 
kiloparsecs from the nucleus for the outflowing material. Other studies have shown that 
large distances such as this appear in several FeLoBAL objects (e.g., SDSS J0838-I-2955; 
Moe et al. 2009). Like many FeLoBALs, QSO0318—0600 is an extremely reddened quasar 
(see discussion in Dunn et al. 2010), where the impact of the dust on the full spectral 
energy distribution (SED) is signihcant. Thus, should the dust source be located between 
the outflow and the nucleus (similar to the torus scenario) then the ionizing flux reaching 
the outflowing gas would be severely diminished. This uncertainty in dust location affects 
the distance determination and in turn the determined energy output by up to a factor of 4 


We take the first steps to ascertain the radial location of dust in a sample of FeLoBALs, 
which has implications for the generation of low ionization species. We begin by discussing 
































5 


the method for spectral analysis in §2 and determine the respective sample for this method 
in §3. In §4, we determine and compare the continunm and emission line reddenings for the 
objects, which provides the radial locations of the dust sources. Finally, we summarize the 
findings and discuss the implications in §5. 


2. Diagnostic Tool to Determine Radial Location 


As mentioned earlier, the sources of dust in LoBAL and FeLoBAL quasars could exist 
in two different radial distance regimes, either interior or exterior to the outflowing gas 
with respect to the nucleus. To explain the commonly observed reddening in FeLoBALs, 
we assert logically plausible sources for the dust. For example, given the unih ed model of 


active galactic nuclei (AGN, 


Antonucci &: Miller 


1985 


Urrv fc Padovani 


19951) . a grazing 


sightline past the edge of the torus in the vicinity of the AGN provides the observed 
reddening, which we refer to as scenario 1. The dust torus would provide the simplest 
simultaneous explanation for the reddening and the necessary decrease in ionizing flux 
required to generate low ionization species in these types of outflows. Other structures do 
exist such as nuclear dust spirals observed in spatially resolved low-z quasars and Seyfert 
galaxies that could potentially provide a reddening near the nucleus. In scenario 2, a source 
not necessarily associated with the AGN such as a galactic dust lane at kpc scale distances 
is the cause of the reddening. Without spatially resolved imaging of the host galaxy we 
cannot distinguish which specific structure is present, but we are able to determine whether 
the dust source lies close to the nucleus or at galactic distances. 


Observationally, the torus in scenario 1 would redden both the continuum source and 
the broad emission lines (BELs) that arise from the broad line region (BLR) close to the 
accretion disk-SMBH system. By comparison, the narrow emission lines (NELs), which are 
generated in the narrow line region (NLR) at a much larger distance (100s of parsecs in 









6 


Reynolds et al 


typical quasars; lPetersorull997l). would rem ain unextincted much like the work on Seyfert 


■8 and 1.9 galaxies by I Trippe et al.lboiol) and on the Seyfert 1 galaxy MCG—6 — 30 —15 by 


19971)). In the second scenario, a large structure at a galactic scale would 


redden the continuum, BLR, and NLR by approximately similar amounts. Furthermore, 
a structure at kpc scale distances would be statistically unlikely to cover only one of the 
two AGN components. Thus, direct comparison of the continuum and NLR reddening 
demonstrates whether the reddening occurs interior (scenario 1) or exterior (scenario 2) 
to the NLR. A similarly reddened NLR compared to the continuum (i.e., E(B—V)ArLK ~ 
E(B—V)cont) suggests that the dust lies exterior to the NLR and cannot be explained by a 
dust torus. A reddening value of the continuum significantly larger than that of the NLR 
implies that the dust is in the vicinity of the AGN and is likely associated with the dust 
torus. 

It is important to note that the BELs should be as reddened as the continuum in 
either scenario, as the BLR lies interior to the torus. Ideally we would use the emission 
ratios of the Balmer BELs in conjunction with the NELs to determine if scenario 1 applies. 


recombination in some objects (e.g. 


Kwan & Krolik 

1981; 

Korista & Goad 

2004) 


primarily stems from high densities associated with the BLR region, which leads to radiative 
transfer and collisional effects that impact the description of recombination present in 
the gas. Because measurements of the Balmer decrement of BELs are potentially poor 
approximations of the extinction, we only compare the extinctions of the continuum source 
and the NEL. 

Unlike the BELs, the NELs are adequately described by case-B recombination as the 
number density is significantly lower (10^ cm~^ < tlm < 10 ^ cm“^), which implies that the 


case-B approximation is valid fjOsterbrock fc Ferland 


200611 . Furthermore, the ratio of Ha 
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to H/3 only spans the range 2.74—2.86 (due to density) at temperatures of approximately 
10,000K, implying that the ratio is relatively insensitive to temperature differences. 


Assuming case-B, the reddening is dehned as: 


E{B-V) 


2.21 

Rho — Rh/3 


log 


2.76 

fHa ’ 
fnp 


( 1 ) 


where Rna and Rnp are the values of the SMC reddening curve aX Ha and Hf3, respectively, 
and fua and fng are the peak fluxes of the emission lines. We use 2.76 as the intrinsic 
Balmer ratio as determined by case-B recombination (Osterbrock & Ferland 2006). While 
Ha and H/5 are the optimal lines for this process, in the majority of our sample we are 
only able to observe H/3 and Hy due to the combination of redshift and spectral range of 
the Sloan Digital Sky Survey (see §3). Thus, we also use the ratio of Hy to H/5 (0.474, 
Osterbrock & Ferland 2006), which likewise spans a relatively narrow range of 0.469—0.476 
for case-B recombination. 


3. Survey 

3.1. Sample Selection 


To employ the method outlined in Section 2, we require a sample of quasars that 
demonstrate FeLoBAL troughs as well as strong narrow line Balmer emission. To obtain the 
largest number of objects, we utilize the Sloan Digital Sky Survey (SPSS ) spectroscopica lly 
obser ved quasars (Richards et. al 2009) through Data Release 7 (DR7, 


Schneider et ah 


20101) . Quasar spectra in the SDSS-III catalog are unusable for 


in that catalog are not accurately flux calibrated flDawson et al. 


lis st udy as the spectra 


20131). We access and 


download the calibrated spectra through the online Catalog Archive Server (CAS). 


We limit our search primarily by redshift. A redshift of z=0.45 is sufficient to redshift 











the Mg II (AA2796,2804) emission line into the spectral range of SDSS (approximately 
3800 to 9200 A) and permit absorption trough detections for Mgll up to approximately 
10,000 km s“^. A redshift of 0.45 places the Balmer Ha (A6563) at ~9500 A and out of 
the spectral range of SDSS, thus precluding simultaneous observation of BAL troughs and 
the Ha emission line. We select an upper limit on the redshift determined by the Balmer 
H/? (A4861) line, which redshifts out of the spectral range at approximately z=0.90. We 
further limit the sample by magnitude as objects fainter than r magnitude of 19 will have 
substantial noise that can potentially obscure weaker narrow Balmer emission lines. The 
resulting sample contains 8,511 quasars. 


Next, we plot and visually inspect these spectra for blueshifted Mgll and Fell BALs. 
To reduce the search time we examine spectra for Fe ll troughs in those objects that contain 
Mgll BALs; previously these troughs have only been observed in objects that show Mgll 
absorption, which is typically easier to identify. We dehne BALs here to be any intrinsic 


absorption troughs in t 


Barlow fc Sargent 


le quasars spectrum (i.e., outflowing gas ejected by the AGN, see 


1997 1. Zhang et al. (2009) suggested that a velocity width of 1600 


km s“^ simultaneously maximizes the frequency of which BALs are intrinsic outflows and 
also maximizes completeness of a BAL sample. To this end, we only include objects with 
contiguous trough widths greater than 2,000 km s“^. A width of 2,000 km s“^ is similar 
to the Absorption Index {AI) method used in Trump et al. (2006) and Hall et al.(2002). 

A width of 2,000 km s“^ is also sufficient to blend the two members of the Mgll doublet 
and facilitate visual inspection. This provides a conservative lower limit on objects with 
intrinsic absorption. Comparing to the Zhang et al. determinations, our sample is complete 
to slightly less than 75%, but over 90% of the objects are bona fide outffows tied to the 
quasar. Quasars in the lower redshift ranges only provide minimal velocity coverage of the 
UV 3 Fell multiplet (~A2600). For the redshift range of 0.45—0.50, we detect only those 
BALs at low velocity (v~6000 km s“^ at z=0.50) and obtain only a lower limit on the 
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number of objects with Fe ll BAL troughs. 


The last criterion we impose on the sample is to select objects with Mgll and Fell 
BALs that show narrow emission lines due to H/3 and/or H 7 (A4341) in conjunction with 
[O III] A5007 emission. We characterize this as a distinct emission feature that has a 
maximum FWHM le ss than 1000 km s“^. This value is consistent with the findings of 


Richards et al 


2003 


for [O ll] A3727 and [O ill] emission lines in reddened quasars. The 


peak intensity must be distinguishable from the BEL (i.e., at least 3 a above the noise level 
at the BEL peak). 


3.2. Survey Objects and Properties 

In our sample, we hnd 70 quasars (approximately 0.8% of the full 8,511 quasar sample) 
that show Mg II troughs matching the criteria established in the previous section. This 
is similar to other spectroscopically selected samples in the optical (e.g., approximately 
1.3% via the AI method in Trump et al. 2006). Of these 70 objects, 15 showed BALs 
attributed to Fell (0.2%). To determine the presence of Fell troughs, we check for radial 
velocity agreement with Mgll as demonstrated in Figure 1 for a high resolution spectrum 
of SDSS J0149—1016. Finally, only 4 of the 15 objects also demonstrate evidence of strong 
narrow emission due to hydrogen, [Olll], and [Oil]. Figure 2 displays the SDSS spectra of 
these 4 objects (and 1 additional object explained below). The Mgll and Fell BAL troughs 
are evident in these spectra. Furthermore, in all 4 cases, we find evidence of troughs due 
to excited, metastable state lines of Fell. We list the SDSS identifiers, r’ magnitudes, 
and redshifts for all 70 of the BAL quasars in Table 1 and indicate which objects in the 
sample show troughs due to Fe ll and strong narrow H/5 emission. Of the 4 FeLoBALs with 
detectable H/5 emission, only SDSS J 0802+5513 also has a detectable H 7 narrow emission 
line. 
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Fig. 1.— Absorption trough identifications for Mgi, Mgii, Mnii and several resonance and excited state 
(noted with an asterisk) Fell lines in the SDSS spectrum of SDSS J0802-I-5513. The primary absorption 
trough, marked with a dashed line, occurs at approximately —1600 km s“^. Each panel is plotted in velocity 
from the respective ions listed in the upper right. Other notable lines are indicated with dotted lines at 
shifts of 1600 km s“^ from their respective rest wavelengths. 
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Fig. 2.— Spectra, plotted in the rest frame, for objects demonstrating BALs due to Mgii, with troughs 
due to Fell, and strong NELs ([Oiii] and H/3). The object SDSS names and SDSS measured redshifts are 
listed in the upper left of each panel. The first 4 spectra were taken by SDSS, while the final spectrum 
(SDSS J0149—1016) was taken with the Magellan MIKE spectrograph. The flux uncertainties are plotted 
beneath each spectrum with black histograms. 
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We also include in our study the object SDSS J0149—1016 (R.A. 01:49:06.736, decl. 
-10:16:49.27). While this object has a redshift of 0.364 and is therefore not included in our 
search statistics, we have a Magellan Inamori Kyocera Echelle (MIKE) spectrum of this 


quasar (also sh own in Figure 2 


(for details see 


Bernstein et al. 


. M 


KE is a high throughput double echelle spectrograph 


20031 ). The spectrum was taken on Sept. 15, 2006 and has 


been reduced with the standard MIKE Interactive Data Language (IDL) reduction tools. 
The continuum flux levels agree with the SDSS spectral observation (Sept. 22, 2001). 
SDSS J0149—1016 was observed with a high resolution spectrograph with the intent of 
determining the outflow’s radial distance. This work is still in progress. In the meantime, 
this spectrum provides an excellent opportunity to compare the reddening determinations 
from the combinations of Ha, H/5, and Hy NELs with the continuum as the MIKE spectrum 
has a larger spectral range (~3600—9300 A) than SDSS that provides coverage of the 
Balmer series as well as the Mgll BAL and Fell absorption troughs (see Figured]). 


To strengthen the continuum reddening determinations in 
spectra with data from the 2 Micron All Sky Survey (2MASS; 


^4, we complement t 


he SDSS 


Skrntskie et al. 


19971). These 


IR data provide a longer baseline for the reddening determinations. We use a matching 
radius of 3" and obtain the J, H, and K magnitudes for each object. We list the 2MASS 
magnitudes for the sample in Table 1. The uncertainties for 2MASS data are typically 
0.10—0.15 magnitudes. 


4. Reddening Determinations 
4.1. Continuum Reddening 

We use the SDSS spectra to determine the continuum reddening rather than SDSS 
photometry because the combination of varying BEL emission line fluxes, and more 
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importantly, the presence of BALs can significantly impact the E(B—V) determinations. 
We begin by correcting for Milky Way extinction (E(B—V)Arty) as determined by the onlin e 


extin ction calculator at the NASA/IPAC Extragalactic Database fjSchlaflv fc Finkbeiner 


201ll) . Next, we correct for the intrinsic reddening of the object (E(B—V)q) in the restframe 


of the quasar using an SMC extinction curve (see Section 1; Richard et al. 2003; Hopkins 
et al. 2004). Because the SMC curve may not provide the most accurate description of the 
reddening at short wavelengths (i.e., in the far UV regime; see Dunn et al. 2010), we £t the 
restframe spectrum in the near UV and visible regimes, which are less dependent on the 
dust grain size distribution. To represent a typical un reddened quasar spectrum, we use the 


SDSS composite spectrum (jVanden Berk et al. 


200lh . With regards to the intrinsic shape 


of the full spectral energy distribution (SED), Vanden Berk et al. (2001) and Richards et 


al. (2003) found that the distribution of quasar SED s’ 
with a relatively small spread {a = 0.30). In addition. 


rapes from SP S S da 


Hopkins et al. 


a is Gaussian 


(2004) showed that 


large reddening (E(B—V) > 0.1) occu red in less than 1% o f thei r sample of SPSS obj ects, in 


agreement with other surveys such as 


Urrutia et al. 


(120091) and 


Glikman et al. 


fc012h .Tlius. 


any large deviation from the composite SDSS spectrum is likely due to dust extinction. 


We de-redden each LoBAL spectrum in the sample by the SMG curve to match the 
continuum levels of the SDSS composite spectrum. We complement the SDSS spectrum 
with 2MASS photometry, when available, which provides data over a larger wavelength 
baseline that are less affected by reddening for the fit. The best £t is determined by 
matching regions of continuum between the two spectra shortward of approximately 5000 
a. Due to the differences in the intrinsic continuum slope, presence of emission lines, and 
various depths and velocities of the BALs, we determine hts via visual inspection. While 
these hts are subject to the possiblity of host galaxy contamination in the quasars’ spectra, 
a lack of strong stellar absorption troughs in the restframes of the quasars suggests that 
the host galaxy contributions are relatively small. We derive limits on uncertainty from 
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under and over compensating for the reddening given the limits of the spectral noise in 
the LoBAL’s spectrum. We also model the general shape of the continuum t o prov ide 


predictions of flux levels in the infrafed range with both a 


Mathews k. Ferlandl (119871) SED 


and the “UV-soft” SED of Dunn et ah (2010) and hnd little difference in the hts between 
the two SEDs across the observed range. We show examples of the FeLoBAL hts in Figure 
Eland list the E(B—V)a^vv values, determined quasar extinction values (E(B—V)q), and 
uncertainties for the LoBALs in Table 2. Regarding SDSS J0149—1016, because we have 
the Magellan MIKE and SDSS spectra we measure the continuum reddening for both. We 
list the SDSS measured value (0.10±0.02) in Table 2 and hnd a consistent value from the 
Magellan spectrum of 0.11±0.03. 


As stated in Section 2, the AGN SED slopes likely have a Gaussian distribution 
centered around the slope of the average spectrum. This implies, especially in a sample 
of non-BAL objects, that a fraction of the objects will have bluer slopes than the SDSS 
composite spectrum. We observe this in a small, randomly selected subset of non-BAL 
quasars using our htting technique as a negative E(B—V)q, which suggests that these 
objects have little or no reddening. Similar to the non-BAL quasars, we only hnd a few 
objects in the LoBAL sample with appreciable negative E(B—V)q values. We summarize 
the relative percentages of LoBALs in several reddening ranges in Table 3. Approximately 
49% of the LoBALs in the sample have an E(B—V)q consistent with zero and 23% of the 
objects have signifcant reddening (E(B—V)q >0.1). This is a substantially larger fraction 
compared to a general sample of non-BAL quasars (1% in the same range; Hopkins et ah 
2004) and is consistent with the SDSS Data Release 5 sample of Mgll LoBALs measured 
by Zhang et ah (2010). 
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Fig. 3. — Spectral correction for both Milky Way and host galaxy reddening for the FeLoBAL objects 
in the SDSS survey. The black histogram is the SDSS spectrum, the red histogram is the SDSS composite 
spectrum, the red curve is the “UV-soft” SED adjusted to the SDSS flux at 5100 A, the green histogram is 
the dereddened spectrum, and the blue asterisks and diamonds are the flux values and uncertainties derived 
from the 2MASS and UKIDSS photometry, respectively. The SDSS composite spectrum likely deviates from 
the SED ( and the 2MASS data) lo ngward of 5000 A due to host galaxy contamination in the composite 


spectrum (|Vanden Berk et al 


2001ri . The derived values of E(B—V) are listed along with the redshift of 


each object in the upper right. Note that the data near 4500 A for SDSS J08 02+5513 are poor matc hes to 


the SDSS spectrum, primarily due to strong Fell emission (similar to EVl of 


Boroson fc Green 


199211 . The 


poor fit to the UKIDSS data for SDSS J1344+3317 is plausibly due to a 10% decrease in flux of the object 
in the time span between the observations (approx. 2 years) given the similar slopes. 
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4.2. NLR Reddening 


To determine the reddening of the hydrogen NELs, we hrst normalize the continuum 
and the BEL (as we do not require any physical information from the BLR, see §2) with a 
spline £t and subtract the fit to isolate the NEL. Using the [O ill] A5007 emission line as 
a template for the NEL, we scale the template to match the normalized hydrogen NELs 
(as demonstrated in Figure 0] for SDSS J0802+5513). We use this template primarily 
to determine the portion of the total emission prohle that is due to the NEL and to 
also account for any asymmetries in the line profile commonly foun d in NELs such as 
multiple peaks (e.g . , SDSS J1344+3317) and asymmetric blue wings flHeckman et al. 


1981 


Greene &: Ho 


2005 


Komossa et ah 


20081) ■ We do not model the lines with Gaussian hts. 


which would require several hts that do not yield any physical information. To determine 
the measurement uncertainties, we add the average hux uncertainty percentage across the 
peak of the line with the continuum determination uncertainty (approximately 3%). These 
are listed in Table 2. 


For SDSS J0149—1016 (also shown in Figure 0]), there are three important features. 
First, we £t Ha, H/3, and Hy due to the broader spectral coverage. Second, this object also 
shows structure in the narrow emission lines, seen in both the [Olll] doublet AA4959, 5007 
and Balmer series. Finally, this object shows blueshifted Balmer absorption troughs, which 
implies a high density for the outhowing gas (see §5 for further discussion of the physical 
implications). In the measurements of the emission features, the absorption troughs are 
important as they directly impact the blue emission wing for fitting. Therefore, our matches 
are primarily compared to the red wings and peaks of the emission lines. 

In the case of SDSS J1632+4204, we are unable to generate an adequate template for 
the NEL. The [O ill] doublet (AA4959, 5007) lies in the telluric absorption corrected region 
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4820 4840 4860 4880 4330 4340 4350 4360 

Rest Wavelength (A) 

Fig. 4.- Left: Plots of the H/3 A4861 and H 7 A4341 narrow emission lines with their fits for SDSS J0802H-5513. The 

black histograms are the hydrogen emission (after subtracting the broad emission lines) and the red histograms are templates 
created from the [O ill] emission line (A5007) scaled to fit the hydrogen lines. The [O ill] lines are boxcar smoothed by 5 pixels, 
which effectively removes problematic pixels due to telluric line subtraction and maintains the general shape of the [O ill] line 
(FWHM = (800 ib 100) km/s). The poor fit to the H 7 red wing is presumably due to contamination by [Oiii] A4363. Thus, we 
fit the blue wing and peak. Residual differences are ploted in the lower sections. We determine uncertainties in fitting by taking 
the difference of matching the peaks to matching the emission line wings, which are comparable to the average of the residual 
differences (better than 20% across several pixels). Right: Plot of the Ha and H/3 emission lines in SDSS J0149—1016. The 
data is rebinned by a factor of two and yields a ratio of 3.23 between Ha and H/3. The data points blueward of the emission 
peak reflect the BALs seen in the Balmer series. The largest differences are due to spurious points in the data, which are from 
imperfect telluric corrections. 


of the spectrum, which provides signihcant errors in an [O ill] template as well as any 
measurement of the H/? emission. There is also no indication of a strong [O ll] A3727 line 
to provide a template. Due to the absence of an uncontaminated line for the template but 
clear presence of H/3 emission, we include this object in the survey statistics, but cannot £t 
the H/3 line. 

As previously identified in Section 3, the only object in the SDSS survey that exhibits 
a measurable narrow emission line from Hy is SDSS J0802+5513. To determine the NLR 
reddening values in the remaining two objects with measurable H/3 emission lines, we use 
the template to determine the upper limit of Hy emission based on the statistical error of 
the data (similar to the method used in Dunn et ah 2010 for determining absorption trough 
limits). We compare a scaled emission line template to the continuum region where the line 
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would be detected. We maximize the scale to match the amount of noise present in the 
spectral region, which provides the largest possible amount of emission present from H7. 
We show this for SDSS J1344+3317 and SDSS J2107—0620 in Figure El Both objects show 
some signs of a weak and noisy H7 line, but we conservatively list these measurements as 
limits. Finally, we list the E(B—V)ArLij determinations and limits for these two objects in 
Table 2 for each object in the sample. 


4 
3 
2 

1 

0 

1 
0 

-1 

4840 4860 4880 4900 4340 4350 4360 4370 4855 4860 4865 4870 43304335434043454350 




Rest Wavelength (A) Rest Wavelength (A) 

Kig. 5. - Left: Similar plot to Figure [4] however for SDSS J1344+3317. The fit to the H 7 emission, is a firm upper 

limit determined by the flux uncertainties. As with Figure [4] the excess emission near 4370 A is likely due to [Oiii] A4363. 
No boxcar smoothing was applied to the [Oiii] template for this object, as the emission lines in this object are free of telluric 
contamination. Right: Plot for SDSS J2107—0620. Unlike with SDSS J1344+3317, we boxcar smooth the oxygen line profile 
by a factor of 3 as the line is much noisier in this object’s spectrum. Also peculiar in this spectrum, the oxygen line appears to 
be shifted by approximately 2 A from the Balmer lines, but does not affect the overall determined ratio. 


While the major focus of this study pertains to FeLoBALs due to the recent distance 
determinations for the outflows of this type and the problems associated with reddening, we 
have also measured the NEL reddening for the regular LoBALs using the techniques listed 
above for the FeLoBALs. We list the values of E(B— V)ArL_R for the LoBALs in Table 2. 
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4.3. Reddening Comparisons 


Many studies of NEL/BEL/continuum reddening in AGN (e.g., Reynolds et al. 1997) 


compare 

However, 


lydrogen column deii sities {Nh) between components of the unified model. 


Maiolino et al. 


(120011) showed that the E(B—V)/A'i^ ratio ranges by two orders 


of magnitude depending on the type of dust in the object. It has also been shown (Dunn 
et al. 2010; Hall et al. 2004) that while the reddening curves in many AGN are best fit 
by an SMG curve, in extreme cases of reddening the continuum is not well fit at short 
wavelengths. This disparity potentially affects any dust column density determinations. 
Thus, we directly compare E(B—V) values between the AGN components. 


Gomparing the values for the FeLoBALs in Table 2, we hnd that three objects show 
clear evidence of reddening in both the continuum and the NELs (SDSS J0802+5513, 
SDSS J1344+3317, and SDSS J2107-0620). While SDSS J0149-1016 shows measurable 
continuum reddening, due to the uncertainties in the line fitting, the ratios and the resulting 
E(B—V)Ar£;L values are consistent with either reddening scenario. As stated previously, the 
NELs in SDSS J1632+4204 cannot be measured, but the object does show an appreciable 
continuum reddening. 


All three FeLoBALs with measurable H/1 NELs listed above, demonstrate a dusty 
source that simultaneously reddens both the NLR and the continuum. This implies that 
the dust must exist radially exterior to the NLR and precludes a torus reddening source 
for these three FeLoBALs. Due to the relatively low number of available objects and the 
combination of low continuum reddening (E(B—V)q ~0.1) with large uncertainties in 
SDSS J0149—1016, we cannot state conclusively if FeLoBALs in general share this physical 
picture. It is notable, though, that all three objects with significant continuum reddening 
appear to have reddening sources farther from the central black hole than the NLR. 

For the regular LoBALs, the majority of the LoBALs have E(B—V)q values consistent 
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with no reddening; we find that the majority of the E(B—V) nlr values are likewise small. 
There are 6 LoBALs of the 70 in the survey with E(B—V)q > 0.1. Of these 6 objects, three 
have measurable H/5 lines: SDSS JlOlO+1843, SDSS J1700+3955, and SDSS J1703+3839. 
We find an E(B—V)q= 0.11±0.03 for SDSS JlOlO+1843. The emission line measurements 
yield a similar value though with a significant uncertainty. SDSS J1700+3955 has moderate 
continuum reddening with a relatively large uncertainty (E(B—V)q= 0.16±0.14). The 
measured NLR reddening is consistent with the continuum within the large uncertainties of 
the two values (E(B—V)q= 0.59±0.30). Finally, SDSS J1703+3839 has a large reddening 
value of E(B—V)q= 0.40±0.03, but unfortunately has no measurable Hy emission line. 

The limit derived from Hy yields an E(B—V)7VLfl >0.2 and is certainly in line with the 
continuum determination. Simlar to the FeLoBALs, the LoBALs in the survey also have 
continuum and NLR reddenings that are similar within the uncertainties, which tentatively 
suggests that the reddening sources in LoBALs exist at larger distances than their NLRs. 

To illustrate the relationship, we plot in Figure |6] the continuum reddening against 
the NLR reddening for both the LoBAL and FeLoBAL quasars in the entire sample. The 
plot shows a significant agreement between the two reddenings determinations for both 
populations of quasars. This nearly 1-to-l relationship supports scenario 2 for the physical 
picture for the quasars in this sample. 


5. Summary and Conclusions 

We have assembled a subsample of 70 LoBAL quasars from a sample of 8,511 quasars 
with redshifts between 0.45>z>0.80 from the SDSS. Of the 70 quasars, 23 have detectable 
narrow emission lines from H^. This sample also contains 15 FeLoBAL quasars, 4 of 
which show signs of narrow Balmer emission lines. Therefore, narrow Balmer emission only 
appears to arise in approximately 27% of FeLoBAL quasars. Both the LoBAL sample and 
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E(B-V)^ 


Fig. 6.— Comparison of the reddening for the quasar’s continuum (horizontal axis) and the reddening of 
the NLR (vertical axis). Diamonds represent the FeLoBALs, while x’s represent LoBALs. Vertical arrows 
illustrate objects with a lower limit determination for the NLR reddening. 

FeLoBAL subsample clearly have a large fraction of objects with signihcant continuum 
reddening (~23% with E(B—V)q > 0.1 for Mgll BALs and 4 Fell objects or ~23% with 
E(B-V)q > 0.1). 

From measurements of the reddening of the continuum and narrow line region, we hnd 
that three FeLoBAL quasars (out of 5 total objects, 4 from our sample and one additional 
quasar with a Magellan MIKE spectrum) have reddening sources that must exist radially 
exterior to the NLR from the AGN and are not directly tied to the dust torus. Due to 
telluric contamination, the H/3 Balmer emission in one object (SDSS J1632+4204) cannot 
be measured. Although the hnal object shows NLR reddening that is consistent with a 
radially exterior source, due to relatively low reddening (E(B—V)q RiO.l) and signihcant 
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uncertainties in the emission line fits, we cannot conclusively state which physical picture 
applies. Also, as strong NELs do not appear to be ubiquitous to FeLoBALs, these objects 
could potentially be inherently different from the general population. 

Given that every measurable object in the FeLoBAL sample has a similar physical 
picture, the Endings hint that the reddening frequently observed in FeLoBALs is due to a 
source situated outside of the outflow from the AGN. In studies of FeLoBALs, (e.g., Moe et 
al. 2009; Dunn et al. 2010; Borguet et ah 2013) the outflow distance determinations have 
all yielded distances of several kpc. Because the objects in our FeLoBAL survey all suggest 
that the dust is located radially exterior to the outflow, the dust would likewise exist at 
kpc scale distances (assuming these objects are similar to previously studied FeLoBALs). 

A simple explanation for a host of the dust would be the outflows themselves. This could 
be further evidence that FeLoBAL quasars are evolutionary states of disrupted systems 
associated with major galaxy mergers as suggested by Glikman et al. (2012) and others. 

In the case of SDSS J0149—1016, should the dust lie at a galactic scale, the implied 
outflow distance is likely significantly smaller than other measured FeLoBALs. This 
distance is implied by the presence of hydrogen Balmer line troughs, which are uncommon 
in BAL quasars. Hydrogen Balmer absorption typically suggests a rather large number 
density (albeit this depends on the ionization state of the gas). Thus, the dust would not 
arise due to the outflow but due to a source beyond the NLR such as a galactic dust lane 
or a nearby dwarf galaxy not associated with the outflow itself in at least this one object. 

We have also examined the continuum and NLR reddening values in the larger sample 
of LoBAL quasars. While many of these objects have low continuum reddening values 
(E(B—V)q <0.1), three objects exhibit relatively large reddening values. The two largest 
reddened objects both have Balmer emission lines that appear to have similarly large 
reddening; this suggests that like the FeLoBALs the dust exists at larger distances than the 
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NLR for LoBAL quasars. 

Due to the rare nature of LoBAL and FeLoBAL quasars and the subsequent small 
sample size (especially with regards to the FeLoBAL quasars), the next step is to increase 
the number of targets. Obtaining spectra for higher redshift FeLoBALs observed by SDSS 
with higher resolution and larger wavelength coverage (similar to SDSS0149—1016) will 
increase the sample size and test these hndings. Infrared spectra will provide coverage of 
the restframe optical regime for the stronger Ha and H/I emission lines and yield a greater 
number of FeLoBALs with narrow Balmer emission suitable for the technique outlined in 
this study. 

LCH acknowledges additional support from the Kavli Foundation, Peking University, 
the Chinese Academy of Sciences, and the Carnegie Institution for Science. 



Table 1. SDSS Quasars with Mgll Outflows and Properties 


SDSS Name 

Plate-MJD-Fiber“ 

r’ Mag 

J Mag'^ 

H Mag'^ 

K Mag'^ 


Fell'^ 

H/3" 

J010352.46+003739.7 

51816-0396-471 

17.42 

16.157 

15.425 

14.074 

0.703 



J024220.10-085332.7 

51910-0456-291 

18.86 

- 

- 

- 

0.799 


y 

J080248.19+551328.8 

53384-1871-440 

18.42 

16.434 

15.923 

14.645 

0.663 

y 

y 

J080611.02+115029.0 

53794-2418-616 

18.88 

- 

- 

- 

0.616 

y 


J080934.64+254837.9 

52670-1205-588 

17.93 

17.029 

16.469 

16.099 

0.545 


y 

J082231.53+231152.0 

53317-1926-546 

17.68 

- 

- 

- 

0.653 


y 

J083525.98+435211.3 

52232-0762-175 

17.55 

16.407 

15.961 

15.131 

0.568 



J085053.12+445122.4 

52605-0897-359 

17.41 

16.222 

15.967 

14.909 

0.542 



J085357.88+463350.6 

52238-0764-248 

18.22 

16.758 

15.892 

15.287 

0.549 



J092525.16+202139.0 

53708-2289-183 

18.51 

17.067 

16.080 

15.288 

0.460 



J093228.56+555344.8 

51991-0556-251 

18.76 

16.611 

15.820 

14.826 

0.876 



J093315.09+351944.2 

52992-1594-348 

18.87 

- 

- 

- 

0.525 


y 

J094225.42+565613.0 

52253-0557-359 

18.52 

- 

- 

- 

0.831 



J094443.13+062507.4 

52710-0993-535 

16.24 

14.831 

14.517 

13.808 

0.695 





Table 1—Continued 


SDSS Name 

Plate-MJD-Fiber“ 

r’ Mag 

J Mag'^ 

H Mag'^ 

K Mag'^ 


Fe ll'^ 

J101038.76+184321.8 

53768-2373-124 

18.97 

16.998 

16.123 

15.433 

0.745 

y 

J101729.52+264146.7 

53765-2350-435 

18.70 

- 

- 

- 

0.526 

y 

J102802.33+592906.7 

52316-0559-582 

18.92 

- 

- 

- 

0.535 


J103036.93+312028.8 

53440-1959-121 

17.69 

16.368 

15.815 

15.059 

0.873 


J103255.37+083503.2 

52734-1240-316 

17.73 

15.750 

15.641 

15.229 

0.891 

y 

J104122.84-005618.4 

51913-0274-060 

18.50 

16.951 

16.176 

15.358 

0.497 

y 

J104210.43+501609.1 

52354-0875-583 

18.07 

16.650 

16.238 

15.600 

0.787 


J104459.60+365605.1 

53463-2090-329 

16.81 

15.537 

15.157 

14.331 

0.701 

y 

J105259.99+065358.0 

52670-1001-080 

18.56 

- 

- 

- 

0.722 


J105856.75+480805.4 

52646-0964-305 

18.00 

16.539 

15.965 

15.313 

0.591 

y 

J111628.00+434505.8 

53061-1364-095 

17.44 

15.887 

15.748 

15.291 

0.801 


J112526.12+002901.3 

51614-0281-427 

18.10 

16.485 

16.300 

15.356 

0.864 

y 

J112621.25+343628.9 

53713-2100-539 

17.86 

16.545 

16.308 

15.594 

0.603 


J112822.42+482310.0 

52642-0966-203 

17.59 

16.546 

16.265 

15.919 

0.543 




Table 1—Continued 


SDSS Name 

Plate-MJD-Fiber“ 

r’ Mag 

J Mag'^ 

H Mag'^ 

K Mag'^ 


Fell'^ 

R/3^ 

J112828.31+011337.9 

51992-0512-123 

18.60 

16.055 

15.386 

14.440 

0.893 

y 


J113807.83+531231.6 

52367-0880-404 

18.67 

16.953 

16.406 

16.107 

0.790 



J114043.62+532438.9 

52734-1015-085 

18.31 

16.942 

16.408 

15.304 

0.530 



J114209.01+070957.7 

53383-1621-306 

18.42 

- 

- 

- 

0.497 



J120751.48+253953.7 

54484-2656-420 

18.18 

16.579 

16.361 

15.802 

0.874 


y 

J121113.38+121937.3 

53149-1612-179 

18.60 

16.764 

15.894 

14.893 

0.464 


y 

J121303.40-014450.9 

52367-0332-579 

18.54 

- 

- 

- 

0.612 


y 

J121442.30+280329.1 

53823-2229-557 

17.38 

15.790 

15.649 

14.594 

0.695 

y 


J123820.19+175039.1 

54234-2599-503 

16.42 

15.326 

14.285 

12.957 

0.453 


y 

J124300.87+153510.6 

53502-1769-584 

18.52 

- 

- 

- 

0.561 

y 


J130741.12+503106.4 

52753-1281-361 

18.25 

- 

- 

- 

0.701 



J131433.19+471457.7 

53062-1461-431 

18.88 

15.518 

15.017 

14.606 

0.869 



J134415.75+331719.1 

53503-2024-346 

18.82 

- 

- 

- 

0.686 

y 

y 

J134651.31+302421.7 

53851-2094-598 

18.44 

_ 

_ 

_ 

0.867 





Table 1 



H Mag^ K Mag' 


z 


Fe ll'^ H/3® 


15.646 

14.993 

0.830 

16.015 

15.176 

0.530 

15.585 

14.845 

0.594 

- 

- 

0.863 

14.903 

14.322 

0.801 

- 

- 

0.587 

16.074 

15.447 

0.717 

- 

- 

0.637 

16.630 

15.491 

0.761 

15.696 

14.859 

0.798 

- 

- 

0.703 

14.866 

13.860 

0.661 

15.251 

14.892 

0.835 


14.301 13.224 0.849 



Table 1—Continued 


SDSS Name 

Plate-MJD-Fiber“ 

r’ Mag 

J Mag'^ 

H Mag'^ 

K Mag'^ 


Fell'^ 

R/3^ 

J154620.98+453916.7 

52782-1333-601 

17.57 

16.274 

16.134 

15.594 

0.459 


y 

J160143.75+150237.7 

54568-2524-356 

17.26 

15.772 

15.638 

14.587 

0.650 

y 

y 

J160234.88+160041.1 

53555-2197-225 

18.18 

- 

- 

- 

0.719 



J160656.95+133931.2 

54569-2527-346 

18.67 

- 

- 

- 

0.452 


y 

J161637.15+390356.8 

52759-1336-088 

17.90 

16.970 

16.333 

15.344 

0.810 



J162435.28+090731.6 

54589-2532-378 

18.06 

16.547 

15.936 

14.949 

0.652 


y 

J163255.46+420407.8 

52379-0816-569 

18.68 

- 

- 

- 

0.726 

y 

y 

J163656.84+364340.4 

52782-1174-337 

18.91 

17.013 

16.251 

15.719 

0.850 



J164447.19+311437.2 

52781-1340-002 

17.90 

16.369 

16.115 

15.281 

0.690 



J170010.82+395545.8 

52079-0633-482 

18.99 

17.240 

16.463 

15.641 

0.577 


y 

J170341.82+383944.7 

52071-0632-632 

18.72 

16.795 

15.999 

15.426 

0.554 


y 

J204333.20-001104.2 

52435-0981-044 

17.94 

16.793 

15.769 

15.223 

0.545 



J210757.67-062010.6 

52174-0637-610 

17.22 

15.732 

15.082 

14.145 

0.642 

y 

y 

J220931.92+125814.5 

52519-0735-501 

18.52 

- 

- 

- 

0.813 


y 



Table 1—Continued 


SDSS Name Plate-MJD-Fiber“ r’ Mag J Mag'^ H Mag'^ K Mag'^ Fell^ H/3® 


^The corresponding SDSS plate number, modified Julian date, and fiber number for the object, 
'’redshift of the object. 

'^IR magnitudes from 2MASS. 

'^Absorption troughs from Fe ll detected at similar velocities to Mg ll. 


BALs with narrow H/3 emission lines. 



Table 2. LoBAL Reddening Determinations 


SDSS Name 

E(B—V)mw 

E(B-V)g 

E(B—V)Ar2,R 

FeLoBALs 

0149-1016“ 

0.038 

0.11+0.03 

0.15+0.15 

0802+5513 

0.042 

0.35+0.05 

0.62+0.26 

0806+1150 

0.022 

0.10+0.11 

— 

1032+0835 

0.021 

0.08+0.07 

— 

1044+3656 

0.013 

0.05+-0.04 

— 

1125+0029 

0.031 

0.11+0.09 

— 

1128+0113 

0.029 

0.10+0.10 

— 

1214+2803 

0.021 

0.08+0.07 

— 

1243+1535 

0.032 

0.00+0.10 

— 

1344+3317 

0.014 

0.20+0.10 

>0.23 

1442+5336 

0.009 

-0.05+0.05 

— 

1448+4043 

0.011 

-0.02+0.09 

— 

1513+2002 

0.032 

0.11+0.05 

— 

1601+1502 

0.042 

0.06+0.04 

— 

1632+4204 

0.025 

0.11+0.09 

_b 

2107-0620 

0.058 

0.14+0.05 

>0.18 

LoBALs 

0103+0037 

0.029 

0.13+0.06 

— 







Table 2—Continued 


SDSS Name 

E(B—V)mw 

■ E(B-V)q 

E(B—V)p^lr 

0242-0853 

0.022 

0.03+0.07 

0.06+0.28 

0809+2548 

0.032 

-0.02+0.09 

> -0.05 

0822+2311 

0.034 

0.00+0.03 

> -0.07 

0835+4352 

0.028 

0.00+0.04 

— 

0850+4451 

0.024 

-0.03+0.07 

— 

0853+4633 

0.021 

0.00+0.03 

— 

0925+2021 

0.038 

-0.05+0.15 

— 

0932+5553 

0.025 

0.15+0.05 

— 

0933+3519 

0.012 

0.00+0.10 

-0.06+0.18 

0942+5656 

0.017 

0.06+0.04 

— 

0944+0625 

0.026 

0.03+0.03 

— 

1010+1843 

0.027 

0.11+0.03 

0.09+0.28 

1017+2641 

0.023 

0.00+0.10 

> 0.03 

1028+5929 

0.007 

0.00+0.10 

— 

1030+3120 

0.018 

0.00+0.05 

— 

1041-0056 

0.046 

-0.10+0.03 

0.14+0.16 

1042+5016 

0.013 

0.03+0.05 

— 

1052+0653 

0.033 

0.00+0.09 

— 

1058+4808 

0.011 

-0.02+0.04 

0.18+0.14 

1116+4345 

0.013 

0.06+0.04 

— 






Table 2—Continued 


SDSS Name 

E(B—V)Miy 

E(B-V)q 

E(B—V)p^lr 

1126+3436 

0.022 

-0.03+0.03 

_ 

1128+4823 

0.017 

0.05+0.04 

— 

1138+5312 

0.007 

0.11+0.09 

— 

1140+5324 

0.009 

-0.05+0.15 

— 

1142+0709 

0.034 

0.00+0.10 

— 

1207+2539 

0.018 

0.05+0.05 

0.03+0.19 

1211+1219 

0.025 

0.06+0.14 

0.13+0.26 

1213-0144 

0.020 

0.03+0.04 

0.27+0.21 

1238+1750 

0.021 

-0.04+0.04 

-0.05+0.12 

1307+5031 

0.013 

0.08+0.07 

— 

1314+4714 

0.009 

0.03+0.07 

— 

1346+3024 

0.017 

0.00+0.03 

— 

1408+3054 

0.009 

0.05+0.05 

— 

1426+0325 

0.031 

-0.02+0.07 

— 

1429+5238 

0.012 

0.02+0.08 

— 

1452+2500 

0.030 

-0.03+0.13 

— 

1457+4521 

0.015 

0.02+0.05 

— 

1457+5234 

0.016 

-0.10+0.11 

— 

1458+4330 

0.014 

0.08+0.12 

— 

1508+3404 

0.014 

0.06+0.04 

-0.19+0.17 






Table 2—Continued 


SDSS Name 

E(B—V)mw 

■ E(B-V)g 

E(B—V)Ar2,ij 

1523+3914 

0.018 

0.04+0.03 

— 

1532+0613 

0.042 

0.02+0.03 

— 

1543+1624 

0.028 

0.06+0.04 

— 

1546+4539 

0.013 

0.00+0.10 

-0.05+0.24 

1602+1600 

0.028 

0.02+0.05 

— 

1606+1339 

0.036 

0.03+0.12 

0.02+0.12 

1616+3903 

0.007 

0.00+0.10 

— 

1624+0907 

0.057 

0.00+0.10 

0.02+0.90 

1636+3643 

0.012 

0.10+0.10 

— 

1644+3114 

0.024 

0.23+0.07 

— 

1700+3955 

0.020 

0.16+0.14 

0.59+0.30^^ 

1703+3839 

0.036 

0.40+0.40 

>0.2 

2043-0011 

0.056 

-0.02+0.08 

— 

2209+1258 

0.072 

0.05+0.05 

>-0.09 


^Object also has Ha. This ratio is between H/3 to 
Ha. 


’^Object has no uncontaminated O ill emission tern- 
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plate available. 

'^Spectrum too noisy for acurate measurement. 
'^Used H/il for line template. 
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Table 3. Reddening Statistics 


E(B-V)q Range 

LoBAL % 

E(B-V) < -0.10 

0.0 

-0.10 > E(B-V)< -0.05 

7.1 

-0.05 < E(B-V) < 0.05 

48.6 

0.05 < E(B-V) < 0.09 

25.7 

E(B-V) > 0.10 

22.9 
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